BCS pairing in fully repulsive fermion mixtures 
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We consider a mixture of two neutral cold Fermi gases with repulsive interactions. We show that 
in some region of the parameter space of the system the effective attraction between fermions of the 
same type can appear due to the exchange of collective excitations. This leads to the formation of 
BCS pairing in the case where bare inter-atomic interactions are repulsive. 



I. INTRODUCTION 

The physics of ultra cold atoms has attracted much at- 
tention in recent decades. One of the reasons is the wide 
variety of quantum collective phenomena which can be 
experimentally observed in these systems. Besides, cold 
atoms can be used as an outstanding playground for in- 
vestigation of condensed matter analogies of phenomena 
belonging to the domain of cosmology, allowing to sim- 
ulate the behavior of such objects as neutron stars and 
black holes f^t— 

One of the main challenges in atomic physics was the 
observation of a Bosc-Einstein condensate (BEC) - a 
macroscopically occupied state formed by bosonic atoms 
cooled down beyond the critical temperature which 
was achieved for the first time in 1995£~— This discov- 
ery stimulated activity in the field of cold bosons and led 
to the experimental investigation of such phenomena as: 
quantized vortices and vortex lattices^ solitonsjiS opti- 
cal trapping^ the Josephson effect and self trapping^ 
the BKT transition^ and others. 

On the other hand, it was realized that cold fermions 
can also reveal interesting physics. Major attention 
was attracted to the BEC-BCS crossover driven by the 
phenomenon known as Feshbach resonance^ It origi- 
nates from the tuning of atomic multiplets with mag- 
netic field and corresponding change of sign of the scat- 
tering length for atoms at the critical magnetic field Bo J£ 
Below the Feshbach resonance, the scattering length is 
positive but individual fermionic atoms are bound into 
molecules, which have bosonic properties and at low tem- 
perature form BEC On the other hand, for B > Bo the 
molecules unbound, but the scattering length becomes 
negative. In this regime the analog of a Bardeen-Coopcr- 
Schricffcr (BCS) many-body state is formedJ^ The latter 
was first proposed for the description of the phenomenon 
of superconductivity which occurs in certain metals and 
alloys and is characterized by formation of loosely bound 
Cooper pairs. Finally, in the vicinity of critical magnetic 
field, Bo, the scattering length diverges which leads to 
strong correlations in the system and onset of the uni- 
versal behavior J£ 

It was recently supposed that the physics of the BEC- 
BCS crossover can be important for understanding of 
mechanisms of unconventional high-tcmpcraturc super- 
conductivity (HTS). The phenomenon was detected ex- 



perimentally for a copper compound^ but the corre- 
sponding theory is still lacking] 20 ' 21 This stimulated pro- 
posals of novel methods of BCS state formation, e.g., 
for the hybrid electron-polariton system^ where the ef- 
fective attractive interaction emerges due to electron- 
exciton interactions. 

In this article we present the idea of a new possible 
mechanism of BCS pairing in a cold fermion system. 
We study the two-component Fermi mixture of neutral 
atoms where all bare particle-particle interactions are re- 
pulsive. However, the total effective inter-particle inter- 
action should account for the possibility of virtual col- 
lective excitations analogous to plasmons in solid state 
systems. In a gas of identical fermions this leads to the 
screening of the inter-particle repulsion. Here we show 
that in a binary mixture similar corrections can lead to 
the onset of an effective attraction and subsequent BCS 
state formation. We analyze the required parameters for 
observation of this phenomenon and present the corre- 
sponding phase diagrams. 



II. THE MODEL 

We consider the ground state of binary fermionic mix- 
ture at ultra-low temperature. The particle-particle in- 
teractions are treated using the s-scattering approxima- 
tion, which is applicable if particles are electrically neu- 
tral and do not have a dipole moment. In this case the 
bare interactions in the system are fully characterized by 
three scattering lengths an, 022 and ai2, corresponding 
to: the interactions between two atoms of type 1; two 
atoms of type 2; and between one atom of type 1 and 
one atom of type 2, respectively. 

The matrix elements of the bare interaction between 
particles of the same kind are given byii 



Vn 



47rai 



mi 2 



and for interaction between different particles 

27T012 
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where mi ^ are masses of the particles and (i = 
fniTO2/(?7ii + denotes the reduced mass. Here we 
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FIG. 1: (Color online). Sketch of the system representing a 
two component Fermi mixture with atoms of the first (red cir- 
cles) and second (green circles) kind. The red arrows schemat- 
ically depict the effective interaction between particles V^V ■> 
which is mediated by the density modulations of the Fermi 
gas consisting of particles of type 2. 



assume all scattering lengths to be positive which means 
that all interactions are repulsive. 

Our goal is to calculate the effective interactions be- 
tween atoms of the same type accounting for the possi- 
bility of creation of virtual collective excitations in the 
ground state of the system. The convenient way to do 
this is to use a language of Fcynman diagrams. While the 
general task of accounting for all possible diagrams can 
not be achieved, there are several approximations which 
allow to explain certain aspects of Fermi gas behavior in 
particular limiting cases. Going beyond the mean-field 
theory, the diagrammatic series, which includes interpar- 
ticle interaction, can be often written using the so called 
random phase approximation (RPA) i 23 i 24 The important 
consequence of the RPA is a correct description of the 
bulk plasmons - collective oscillations of the density. In 
the context of the electronic plasma in metals, the ac- 
count of the RPA corrections leads to the screening of 
the Coulomb interaction, which can be expressed as 



Veff = 



_ Vo(q) 
1-Vb(«)n(g,w) £{q,u) : 



V (q) 



(3) 



where Vo(q) is a bare Coulomb interaction and H(q, w) de- 
notes a polarization operator. e(q, u>) = 1 — Vo(q)TL{q, oS) 
is the dielectric function, which is responsible for charge 
screening and in general is frequency dependent. 

In the case of an electron gas in metals, the RPA is 
applicable in the regime of high densities. The condition 
is that the kinetic energy of Fermi gas should overcome 
its potential energy. The kinetic energy of the three di- 
mensional Fermi gas scales as Ekin ~ n 5 / 3 , where n is 
the concentration of the gas. The potential energy of the 
Coulomb interaction scales as Ep l ot <~ n 4 / 3 . One sees, 
that Ekin is a faster growing function of n, which makes 
theoretical description of the dense electron plasma easier 
than description of the diluted electron gas, for which the 
account of correlation corrections becomes a very tricky 
task. 

However, in the domain of the cold fermion atoms 
the situation becomes qualitatively different. Differently 
from the Coulomb case, the interaction between neu- 
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FIG. 2: Diagrammatic representation of the effective inter- 
action between the atoms of the first kind mediated by the 
collective excitations in a two-component Fermi gas. 



tral atoms is short range, and the potential energy of 
the atomic gas scales as the square of its concentration 
Epot = VL 3 n 2 , with V = Aira/m being the interaction 
constant, L 3 denotes volume and m is the mass of a par- 
ticle. As function of n, it grows faster than the kinetic 
energy, which means that the RPA description in this 
case should work better in the low density limit, the con- 
dition which is usually satisfied in the experiments. We 
will therefore use the RPA approximation in our further 
consideration. 

In the case of multi-fermion mixtures the effective in- 
teraction can be written using the generalized dielectric 
functior* 2 ^^ 
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v (i- v nr 



(4) 



where by bold characters Vo and II we denote 2x2 ma- 
trices of the bare interactions and polarization operators 
describing the system, 



Vn = 



n 



Vn Vi 2 

V12 V22 



IIi 

n 2 



(5) 



(6) 



Subscripts 1 and 2 correspond to the fermion flavor in the 
system. Expression Q can be obtained by performing 
summation of all RPA type diagrams shown in Fig. [2] 

The polarization operator for the Fermi gas can be cal- 
culated using the standard rules of evaluation the bub- 
ble diagrams and, for zero temperature, for the 3D case, 
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FIG. 3: (Color online). The model effective averaged poten- 
tial J7ii(cli) plotted for two different ratios r\ (r\ = 3.5 for red 
lines and r\ — 5 for the blue ones). In the inset we plot the gap 
as a function of frequency for r\ — 5. The model parameters 
used for calculations are an = 0,21 = 200ao where ao denotes 
the Bohr radius. 



read 



,28 



U{q,u) = -d(E F )[l + — J (7) 



where the function f(x,Xo) is defined as 



f(x,x )=[l-(xo/x-x) In 



x + x 



X — X 2 + Xq 



(8) 



where we used dimensionless variables x = q/2kp and 
xq = u/AEp, and d(E) denotes the density of states. One 
should note that the polarization operator Q changes 
the sign at a characteristic frequency u c dependent on q. 

Matrix equation @] allows to obtain the effective inter- 
action between particles of the first kind: 



H2V22) + n 2 y? a 



(1 - iiiVii)(i - u 2 v 22 ) - nin 2 v^ 



(9) 



As expected, V^i (<Z) w ) becomes a function of momen- 
tum and frequency. The first term in the numerator 
corresponds to the screened direct interaction between 
particles of type 1, while the second term describes the 
interactions mediated by collective excitation in the sys- 
tem of particles of type 2 proportional to V^ 2 . Tuning 
the interaction between particles of different kind, i. e., 
making V± 2 ^> Vu, one can make the second term dom- 
inant. Therefore, in the low frequency region where 1I2 
is negative it is possible to reach the regime where the 
total effective interaction becomes attractive. In 

this situation one can expect that the atoms of the type 
1 will form the BCS state, analogous to those formed by 
the electrons in a superconductor. The role played by 
the phonons in a superconductor in our case is played by 
collective excitations of the system 2 analogous to bulk 
plasmons in a metal. The properties of the BCS state, 
such as the value of the gap and critical temperature, 
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FIG. 4: (Color online), (a) BCS phase diagram as a function 
of interaction ratio r\ = an I aw and the critical temperature 
T c . (b) The temperature dependence of the gap A plotted for 
77 = 5 which allows to find the critical temperature of BCS to 
normal phase transition. 



crucially depend on the relative strengths of the interac- 
tion between fermions of type 1 and fcrmions of types 1 
and 2, which can be described by the parameter 



V 



Pig 
an 



(10) 



If two fermions interact attractively, even weakly, they 
can form a Cooper pair. The largest contribution to the 
process is given by the fermions close to the Fermi sur- 
face. Therefore, in the 3D case one can compute an effec- 
tive dimensionless interaction Un(u>) between the atoms 
of type 1 by averaging the potential of interaction V^/ 
over the 2D Fermi sphere, analogous to Rcf. [2^] where 
the case of 2D polariton-mediated superconductivity was 
considered. One has: 



2tt 



V^ ff (q,uj)k%sm6de/M, (11) 



where q — y/2k F (l + cos(0)) is an exchanged momentum 
between two fermions lying on a Fermi surface. The nor- 
malization factor J\f — d(Ep) = rnkp/Trh 2 corresponding 



4 



to the density of states evaluated at the Fermi energy is 
added to make the effective potential dimensionless. 

The shape of Un(ui) is shown in Fig. |3]for two different 
values of 77. One can note the following. First, the effec- 
tive interaction potential is frequency-dependent, which 
means that retardation effects play an important role. 
Second, it crucially depends on 77. For small 77, the di- 
rect repulsion between atoms of type 1 dominates, and 
Uh(lo) is always positive. The increase of 77 leads to 
increase of the role played by interactions mediated by 
collective excitations, and a region of negative values of 
Un(ui) develops at small frequencies. This corresponds 
to the onset of effective attraction at large retardation 
times. The situation is fully analogous to those in super- 
conducting materials, where retarded phonon-mediated 
interactions dominate over direct Coulomb repulsion^ 
Naturally, the attraction becomes stronger with increase 
of 77. 
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behavior of the gap function is shown in the inset of Fig|3l 
It is if fully analogous to those observed in superconduct- 
ing solid state systems. The gap is maximal for lo = 0, 
changes sign at some characteristic frequency and then 
decays to zero when lo — > 00. 

Figure Eta) shows the dependence of the critical tem- 
perature T c on the ratio 77 found by solving Eq. (|12p . 
T c was determined as the value of the temperature above 
which the gap equation has only a trivial solution, A = 0, 
in the whole region of lo. The concentrations of the 
fermions of types 1 and 2 were chosen as 774 = 712 = 10 13 
cm -3 , and atoms of type 2 were taken to be four times 
heavier than atoms of type 1. There exists a critical value 
of the parameter 77 below which BCS pairing is absent 
even at T = (t7 c « 3.15 for the values of the parameters 
we consider). The increase of 77 above its critical value 
naturally leads to the increase of the critical temperature. 

For a given value of 77 the maximum value of the gap 
decreases as a function of the temperature and vanishes 
as T — > T Cl as shown in Fig. HJb) for ?; = 5. 

Finally, in Fig. [5] we present the phase diagram in 
the axes of density of the Fermi gas (n = n\ = 712) and 
the interaction ratio 77. The concentration of the gas 
is a relevant parameter since it enters the polarization 
operator and defines the cut-off energy tno c for which the 
effective interaction potential changes from attractive to 
repulsive. For the BCS phase we use color to plot the 
related critical temperature. As it can be expected, the 
critical temperature of the transition to the BCS state 
increases with the increase of the concentration and the 
ratio 77. 



FIG. 5: (Color online). Phase diagram showing the 
BCS/normal phase transition as a function of the concentra- 
tion of particles of type 2, denoted as n, and the interaction 
ratio 77. The colors indicate the critical temperature in the 
gapped state. 

To find the gap in the spectrum connected to emer- 
gence of the BCS state and Cooper pair formation, we 
solve the gap equation in integral for m 22 ' 30 



A(Lo,T) = -h f 

J — t 



Uh(lo -lo')A(lo',T) ta,nh(E/2k B T) 
2E 



die'. 



^ (12) 

where: E = y/(fiw') 2 + E (k); E = fi 2 k 2 /2mi is the 
kinetic energy of atoms of type 1; and T denotes the 
temperature. 



III. RESULTS AND DISCUSSION 

The gap equation is a non-linear integral equation that 
can be solved numerically by an iterative procedure. The 



IV. CONCLUSIONS 



We studied the influence of virtual collective excita- 
tions represented by diagrams of the RPA type on phys- 
ical properties of a two-component mixture of neutral 
fermions with repulsive interactions. We have shown 
that, in the regime where the interaction between par- 
ticles of the same type is much weaker than interaction 
between particles of different types, the total effective in- 
teraction can become attractive. This can lead to BCS 
pairing in the system. The critical temperature of the 
transition to the BCS state is determined by the concen- 
tration and the value of the parameter 77 which charac- 
terizes the ratio of the interaction between particles of 
the same type and particles of different types. 
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